Half-Metallic Silicon Nano wires: Multiple Surface Dangling Bonds and Nonmagnetic Doping 
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By means of first-principles density functional theory calculations, we find that hydrogen-pas sivated ultrathin 
silicon nanowires (SiNWs) along [100] direction with symmetrical multiple surface dangling bonds (SDBs) 
and boron doping can have a half-metallic ground state with 100% spin polarization, where the half-metallicity 
is shown quite robust against external electric fields. Under the circumstances with various SDBs, the H- 
passivated SiNWs can also be ferromagnetic or antiferromagnetic semiconductors. The present study not only 
offers a possible route to engineer half-metallic SiNWs without containing magnetic atoms but also sheds light 
on manipulating spin-dependent properties of nanowires through surface passivation. 

PACS numbers: 71.70.Ej, 73.21.Hb, 73.22.Dj 



It is widely believed that nanowires can be the most promis- 
ing candidates for the basic building blocks of future nano- 
electronics, as they could be essentially useful in engineer- 
ing diverse nanodevices such as field-effect transistors, logic 
gates, DNA sensors, etc. mm 

Thus, the nanowires with 
thin diameters f?, ^ were actively studied both experimen- 
tally and theoretically in the past years. Recently, particular 
attention was paid on the ultrathin silicon nanowires (SiNWs) 
Q B B E Ini [l2il, as they not only have unusual elec- 
tronic properties of scientific interest, but also could be po- 
tentially incorporated into nowadays well-developed Si-based 
electronic technology. 

The properties of ultrathin nanowires can be modified by 
adjusting diameters and orientations, conducting the surface 
passivation, or introducing dopants, etc. For SiNWs, among 
other things, it is known that (i) due to surface reconstruc- 
tion, SiNWs along [100] direction with diameters smaller than 
1.7 nm prefer a square cross section with sharp corners where 
the electrons are quite localized liilEEQ]; (ii) the ul- 
trathin SiNW that is saturated with H atoms on surface can 
show a semiconducting behavior [8], and if it is substitution- 
ally doped with nonmagnetic elements like B or P, the main 
bands of impurities distribute at the top of valence bands llioll ; 
(iii) H-passivated ultrathin SiNWs along [100] direction inter- 
stitially doped by a certain density of transition metals such 
as Co or Cr can have a half-metallic ground state lUsIl : (iv) 
when the H-passivated nanowires or surfaces have an isolated 
or a single row of surface dangling bonds (SDBs), there is a 
band contributed by the SDBs crossing the Fermi level (FL) 
and if the spin polarization is considered, 
the band splits into two discrete bands, respectively, below 
and above the FL [19], leading to two band groups formed in 
the case of multiple SDBs [20]. 

In the present Letter, by means of ah initio calculations 
we have found that, for an ultrathin H-passivated SiNW along 
[100] direction with multiple SDBs at symmetrical positions 
and substitutionally doping B atoms at the center of cross sec- 
tion, a ferromagnetic (FM) configuration is favored as the 
ground state, and such a SiNW is a half-metal with 100% 
spin polarization, namely, it is insulating for one spin direc- 
tion while metallic for the opposite spin direction at the FL. 
To the best of our knowledge, this is the first time to report 



a novel type of half-metallic SiNWs with multiple SDBs and 
nonmagnetic doping, in contrast to the previous study where 
the half-metallic SiNWs were obtained by doping transition 
metals [15]. This poses a possible route to make half-metallic 
SiNWs through manipulating the surface passivation and in- 
troducing nonmagnetic dopants. 

Our calculations were based on the density functional the- 
ory lElh with generalized gradient approximation expressed 
by PBE functional [22], employing norm-conserving pseu- 
dopotentials L23.] and linear combinations of atomic orbitals. 
The SIESTA code Q Q with double-zeta polarized basis 
sets was used to perform spin-polarized calculations. The 
energy cutoff was 180 Ry, and the total energy convergence 
criterion was 10~^ eV. All structures were fully relaxed by 
the maximum force tolerance of 0.015 eVM, and the lattice 
constant along the wire axis was optimized for each structure. 
Among the SiNWs of 25 Si atoms per unit [Fig. 1(a)] with dif- 
ferent SDBs, the separation between two adjacent wires was 
kept as 15.27 (±0.02) A. The Brillouin zone along the axial 
direction was sampled with 15 k-points. 

Figs. 1(a) and (b) depict the ultrathin SiNW along [100] 
direction passivated with H atoms without SDBs (labeled by 
SiNW-H12), where there are 25 Si and 12 H atoms per unit 




FIG. 1: (Color online) The cross section (a) and side view (b) of 
the optimized structures of H-saturated SiNW along [100] direction. 
The atomic positions marked by numbers in (a) are discussed in the 
context, (c) The SiNW with one SDB per unit. Blue (dark gray) and 
orange (light gray) balls represent Si and H atoms, respectively. 



FIG. 2: (Color online) (a) The band structure and (b) the total and 
projected density of states (DOS) of ferromagnetic SiNW-Hll. The 
Si with SDB is marked in pink (gray). The 8 Si atoms in the corner 
including the SDB are enclosed in the inset of (b) and marked blue 
(dark gray). The red (gray) solid and green (light gray) dashed curves 
in band structures represent majority (t) and minority (|) spins, re- 
spectively, throughout the context. The Fermi energy is set to zero. 



cell. In SiNW-H12, the diameter is about 11 A including H, 
and about 8 A excluding H. The positions in y axis of the ad- 
jacent Si atoms at positions 3 and 4 [Fig. 1(a)] differ only by 
0.009 A. Note that the H saturated SiNW-H12 is semicon- 
ducting Isl]. For the SiNW with single SDB per unit [denoted 
by SiNW-Hl 1, Fig. 1(c)], there is one unpaired electron at the 
SDB in each unit, and the total magnetic moment /i is 1.00 
jj^B per unit cell (/i^ is Bohr magneton). In SiNW-Hll, the 
Si atom with SDB [position 3 in Fig. 1(a)] is stretched 0.047 
A outer than the Si at position 4 that is still passivated with H, 
and the energy for generating a single SDB per unit is 3.52 eV. 
We compared different spin configurations by employing a su- 
percell with two unit cells (including two adjacent SDBs), and 
disclosed that the FM configuration is favored as the ground 
state, with the energy of 1 1 .3 meV per unit cell lower than the 
antiferromagnetic (AF) configuration, and 155.6 meV per unit 
lower than the nonmagnetic (NM) configuration. It turns out 
that the SiNW-Hll is a FM semiconductor with an indirect 
gap of 0.59 eV. 

The electronic structure and density of states (DOS) of FM 
SiNW-Hll are presented in Fig. 2. It can be seen that there 
are two split bands just above and below the FL [Fig. 2(a)]. 
The upper band is of minority spins (|), while the lower is 
of majority spins (j), which correspond to the half-filled band 
crossing the FL if the spin polarization is ignored, similar to 
the case at C(OOl) surface |[l9ll . For the bands away from the 
FL, there is merely a small displacement between the spin- 
up and spin-down subbands, which is close to the SiNW-Hl 2. 
The constituents of the two split bands are clearly manifested 
in the projected DOS (PDOS): The Si atom with a SDB con- 
tributes 44.9% to the total DOS (TDOS), which dominates 
overwhelmingly among all atoms, and the 8 Si atoms in the 
corner including the SDB contribute 82.9% to the TDOS, as 
shown in Fig. 2(b). This indicates that the two split bands 
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FIG. 3: (Color online) The band structures of (a) SiNW-HSFM, 
(b) SiNW-H8AF, and (c) SiNW-H8BFM. (d) The DOS of SiNW- 
H8BFM. The structure and spin configuration of SDBs for (e) SiNW- 
H8BFM and (f) SiNW-H8BAF. The dark blue (dark gray), orange 
(light gray, small), wine (solid) balls represent Si, H, B atoms, re- 
spectively; the pink (gray) and green (light gray) atoms denote Si 
atoms with SDB for majority and minority spins respectively. 



are formed primarily by the unpaired electron of SDB and the 
electrons closely around the SDB. 

Now let us look at the SiNW with multiple SDBs. Without 
loss of generality, the case of four SDBs per unit cell located 
symmetrically on each edge of the cross section will be calcu- 
lated, as shown in Figs. 3(a) and (b), labeled by SiNW-H8. In 
comparison to the SiNW-Hl 1 , the FM configuration of SiNW- 
H8 [SiNW-H8FM, Fig. 3(a)] has more split bands, which are 
separated into two groups above and below the FL, and each 
group includes four bands of the same spin species (either spin 
up or down), where there is one band degenerate with another 
in each group. The magnetic moment is 4.00 /j^b per unit 
cell. However, the SiNW-H8FM is merely a metastable state. 
The ground state has an AF configuration [SiNW-H8AF, Fig. 
3(b)], where the SDBs on adjacent edges have opposite spin 
alignments, with an energy lower than SiNW-H8FM by 26.9 
meV per unit cell. The average energy for generating each 
SDB in this case is 3.55 eV. The spin-up and spin-down bands 
of the SiNW-H8AF coincide, and /i is zero, as revealed in Fig. 
3(b). Therefore, SiNW-H8 is an AF semiconductor with an in- 
direct gap of 0.47 eV. In the band structure of SiNW-H8FM 
[Fig. 3(a)], the indirect gap between the bands of majority 
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and minority spins is 0.34 eV, which is so small that it could 
be possible to generate a perfect half-metallic behavior in the 
SiNW-H8 so long as the two conditions are satisfied: (i) The 
FL is slightly shifted upward or downward to intersect the split 
bands of the same spin species; (ii) The FM configuration of 
SDBs should be the ground state. To meet with the conditions, 
it may be convenient to conduct nonmetallic dopings (such as 
B or P). 

By substitutionally doping one B atom at the center of 
the cross section [position 1 in Fig. 1(a)] of the SiNW- 
H8FM (labeled by SiNW-H8BFM), we found that the SiNW- 
H8BFM is the ground state, and the FL moves indeed down- 
ward intersecting the majority spin bands, giving rise to a 
half-metallic behavior, as shown in Fig. 3(c). This is also 
manifested in Fig. 3(d), where the DOS of majority spin 
G^{Ef) at the FL (Ef) is 8.28/eV, while the DOS of mi- 
nority spin Gi{Ef) is zero, leading to the spin polarization 
P = [G^{Ef) - G|(^f)]/[Gt(^f) + Gi{Ef)] = 100%. 
The indirect band gap of minority spin is 0.86 eV. The total 
magnetic moment of SiNW-H8BFM is 3.00/iB per unit cell, 
while the local moment of each Si atom with SDB is 0.51 
(±0.01) /iB, that is 0.1 1 /iB lower than SiNW-H8AF. As there 
are three spin-up bands at different k points in the Brillouin 
zone crossing the FL, it is hard to open a gap due to Peierls in- 
stability [12]. Thus, the half-metalHcity in the SiNW-H8BFM 
is robust. We also examined possible magnetic configurations 
[FM, AF and ferrimagnetic (FI)] in a supercell with either 
one or two unit cells, respectively, together with NM (SiNW- 
H8BNM, spin-unpolarized) calculations, and observed that 
apart from the FM ground state [SiNW-H8BFM, Fig. 3(e)], 
the lowest metastable state is an AF metal [SiNW-H8BAF, 
Fig. 3(f)]. The total energies of SiNW-H8BNM and SiNW- 
H8BAF are 309.2 meV and 21.5 meV per unit cell, respec- 
tively, higher than that of SiNW-H8BFM. Based on the energy 
difference between the SiNW-H8BFM and SiNW-H8BAF, we 
can estimate that at the temperature lower than 25 OK, the half- 
metallic property in the SiNW-H8B can be retained. 

Doping P at the center in a FM configuration of SiNW-H8 
(labeled by SiNW-H8PFM) is similar to the SiNW-H8BFM 
except that the FL moves upward crossing the spin-down 
bands, resulting in the half-metallicity for minority spins. 
However, the ground state for this case is a FI metal, where the 
total magnetic moment fi = O.OGfiB per unit cell, and the total 
energy is slightly lower than SiNW-H8PFM by 3.2 meV per 
unit cell. Compared with the AF ground state of SiNW-H8, 
doping B at the center will decrease both the charge of each 
SDB and the diameter of the SiNW, leading to a FM ground 
state, while doping P brings about inverse changes, making a 
FI ground state. Besides, it is unraveled that in spite of doping 
B, P or not, if there are two SDBs in one corner [at positions 
3 and 8 in Fig. 1(a)] or in one edge (at positions 3 and 7), the 
FM configuration appears not to be the ground state. 

For the SiNW-H8BFM, Fig. 4(a) indicates that the spin 
density (p^ — Pi) distributes identically in each of the four 
corners with SDB, where p^ (pi) is the number density of 
electrons with spin up (down). The dopant B at the center has 
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FIG. 4: (Color online) The contour plot of spin density distribu- 
tion (in unit of IQ-^ electrons/Bohr^) of (a) SiNW-HSBFM and (b) 
SiNW-H8B(2)FM at the cross section through the central atom. The 
gray, white and dark blue (dark gray) balls stand for the positions 
of Si, H and B atoms, (c) The band structure and DOS of SiNW- 
H8B(2)FM. Four Si atoms with SDBs and one B atom are marked 
pink (gray dash, I), dark yellow (gray short dot, II), purple (dark gray 
dot, III), blue (light gray dash-dot, IV) and wine (solid), respectively. 



a homogeneous effect on the symmetrical SDBs, leading to 
the half-metallicity. If doping B takes place at the edge [posi- 
tion 2 in Fig. 1(a)] instead of the center with the FM configu- 
ration [SiNW-H8B(2)FM, Figs. 4(b) and (c)], there are three 
spin-up bands just below the FL, and four spin-down bands 
and one spin-up band just above the FL. The four Si atoms 
with SDBs contribute to the PDOS distinctly, among which 
the one [red Si (I) in Fig. 4(c)] adjacent to the B atom plays 
an overwhelming role. It is also confirmed by the spin den- 
sity distribution in Fig. 4(b), where the spin density of the 1/4 
cross section around Si (I) is nearly zero, totally different from 
those in other three corners. The local moment of Si (I), (II), 
(III) and (IV) is 0.09, 0.58, 0.63 and 0.61 ps, respectively. 
In general, doping B at any position other than the center will 
give a semiconductor. The total energy per unit cell of doping 
B atoms at position 4, 2, 3, 5, 6 and 7 [Fig. 1(a)] is 0.36, 0.52, 
0.70, 0.46, 0.58 and 0.73 eV, respectively, higher than dop- 
ing at position 1 (SiNW-H8BFM), suggesting the center is the 
most preferred doping position for B. 

For the FM configuration, one B dopant per unit will reduce 
one electron below the FL, which, depending on the relative 
positions among dopants and SDBs, could change the band 
structure in two ways. One is that the dopant affects mainly 
one position of SDBs, moving one spin-up band to above the 
FL, as SiNW-H8B(2)FM. The other is that the dopant affects 
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FIG. 5: (Color online) (a) The SiNW-H8BFM under an external elec- 
tric field E. The DOS under (b) Ez = 5.0 V/i along z direction; and 
(c) Ea: = 0.18, (d) 0.19, (e) 1.0 and (f) 1.8 V/i along x direction. 



equally the multiple positions of SDBs, pushing more spin-up 
bands crossing the FL, as the half-metallic SiNW-HSBFM. 
Consequently, in order to produce the half-metallicity, the 
SDBs should be arranged at symmetrical positions and, the 
dopants must be located symmetrically to the SDBs. In addi- 
tion, several cases with larger radius of 61 Si atoms and four 
symmetrical SDBs per unit cell are also checked. The results 
show that substitutionally doping either one B atom at the cen- 
ter or five around the center per unit cell can still generate the 
FM ground state and lead to the half-metallicity with 100% 
spin polarization. 

To examine whether the half-metallic property of the 
SiNW-HSBFM [Fig. 5(a)] is stable under external electric 
fields, we made calculations in presence of longitudinal and 
transverse electric fields, respectively. By applying an electric 
field along the nanowire direction (along z axis), we found 
that the half-metallic behavior of the SiNW-HSBFM still re- 
tains under a field as large as 5.0 V/A [Fig. 5(b)]. When 
an electric field is applied perpendicular to the nanowire axis 
(along X direction), we observed that the half-metallic state is 
kept when the field is less than O.IS V/A [Fig. 5(c)]; when 
the transverse field reaches 0.19 V/A [Fig. 5(d)] or larger, the 
spin polarization decreases dramatically, and the half-metallic 
property is destroyed, which becomes a FM metal. When the 
electric field increases further, say about 1.0 V/A, the half- 
metallic property comes back again [Fig. 5(e)]. When the 
field exceeds l.S V/A [Fig. 5(f)], it becomes a NM metal. 

In summary, a novel type of half-metallic SiNWs along 
[100] direction by introducing multiple SDBs and boron dop- 



ing is predicted by means of the spin-dependent ab ini- 
tio calculations. The obtained results show that the half- 
metallicity in such SiNWs is quite robust against external 
electric fields. Under the circumstances with different SDBs, 
the H-passivated SiNWs can also be FM or AF semiconduc- 
tors. The present findings might be applicable in nanospin- 
tronics and nanomagnetism. 
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